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, ABSTRACT 

■ We report the detection of eclipses in the close white-dwarf- M-dwarf binary star 

' RX J2130. 6+4710. We present lightcurves in the B, V and I bands and fast photometry 

obtained with the three channel CCD photometer Ultracam of the eclipse in the it', 
' g' and r' bands. The depth of the eclipse varies from 3.0 magnitudes in the u' band 

to less than 0.1 magnitudes in the I band. The times of mid-eclipse are given by the 
ephemeris BJD (Mid-eclipse) = 2 452 785.681876(2) + 0.521 035 625(3)S, where figures 
I 1 in parentheses denote uncertainties in the final digit. We present medium resolution 

spectroscopy from which we have measured the spectroscopic orbits of the M-dwarf 
and white dwarf. We estimate that the spectral type of the M-dwarf is M3.5Ve or 



O 
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M4Ve, although the data on which this is based is not ideal for spectral classification. 
We have compared the spectra of the white dwarf to synthetic spectra from pure 
hydrogen model atmospheres to estimate that the effective temperature of the white 
dwarf is T cfi = 18 000 ± 1000 K. We have used the width of the primary eclipse and 
duration of totality measured precisely from the Ultracam u' data combined with 
the amplitude of the ellipsoidal effect in the I band and the semi-amplitudes of the 
spectroscopic orbits to derive masses and radii for the M-dwarf and white dwarf. The 
M-dwarf has a mass of 0.555 ± 0.023 M© and a radius of 0.534 ± O.O53R,0, which is a 
typical radius for stars of this mass. The mass of the white dwarf is 0.554 ± 0.017M Q 
and its radius is 0.0137 ± 0.0014R Q , which is the radius expected for a carbon-oxygen 
white dwarf of this mass and effective temperature. The lightcurves are affected by 
frequent flares from the M-dwarf and the associated dark spots on its surface can be 
detected from the distortions to the lightcurves and radial velocities. RXJ2130. 6+4710 
is a rare example of a pre-cataclysmic variable star which will start mass transfer at 
a period above the period gap for cataclysmic variables. 

Key words: stars: binaries: close - stars: white dwarfs - stars: fundamental param- 
eters - stars: binaries: eclipsing - stars: individual: RXJ2130. 6+4710 



1 INTRODUCTION 



* Based on observations made with William Herschel, Isaac New- 
ton and Jacobus Kapteyn telescopes operated on the island of La 
Palma by the Royal Greenwich Observatory in the Spanish Obser- 
vatorio del Roque de los Muchachos of the Instituto de Astrofisica 
de Canarias. 



RX J2130. 6+4710 is a soft X-ray source which was de- 
tected in a survey of the galactic plane by the ROSAT 
satellite (Triimper 1983, Motch et al. 1991). Motch et al. 
(1997) used optical spectroscopy of stars near the posi- 
tion of these X-ray sources to identify and classify opti- 
cal counterparts in the region of the Cygnus constellation. 
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RX J2130. 6+4710 was identified with a faint star which they 
labelled RXJ2130. 3+4709 about 12arcsec SSE of the bright 
GO star HD 204906 (=SAO 50947, V=8.45). The spectrum 
of RX J2130. 6+4710 shows the broad Balmer lines of a typ- 
ical DA white dwarf together with narrow emission lines in 
the cores in the Balmer lines and the Ca II H&K lines. Also 
visible are molecular bands of TiO from the Me-type com- 
panion star which is also the source of the narrow emission 
lines. They classified RX J2130. 6+4710 as a WD+Me close 
binary star. Motch et al. report that the radial velocities of 
broad Balmer lines and emission lines move in opposite di- 
rections with an orbital period of approximately 12 hours. 
The white dwarf is too cool to contribute to the X-ray flux 
(1.6 x 10 -2 ctss -1 ) which probably arises from the Me star 
companion. The distance from the optical counterpart to 
catalogued position of the X-ray source is almost exactly 
the same as the 90% confidence radius of 26.1 arcsec in the 
position of the X-ray source. 

It is likely that RX J2130. 6+4710 is the same source 
as the 2E 2128.4+4657, a soft X-ray source detected by 
the Einstein satellite which was identified with HD 204906 
by Thompson, Shelton & Arning (1998). The position of 
2E 2128.4+4657 given in the Einstein Observatory catalog 
(Harris et al. 1994) differs by 23 arcsec from the position 
for the optical counterpart of RX J2130. 6+4710 given by 
Motch et al., which is comfortably within the positional un- 
certainty of 42 arcsec for the X-ray source. We note that 
RX J2130. 6+4710 appears in the "Catalogue of cataclysmic 
binaries, low-mass X-ray binaries and related objects" (Rit- 
ter & Kolb 1998) as "J2130+4710" . The corresponding 
name in the SIMBAD database (Wenger et al. 2000) is 
RX J2130. 6+4710, which is a more accurate reflection of 
the position of the optical counterpart (a = 21 h 30 m 18.6 s , 
8 = +47°10'08", J2000.0) and is the name we use through- 
out this paper. 

White dwarf stars are the remnants of stars less mas- 
sive than about 5-8M© which have been through the red 
giant phase of their evolution. The separation of a WD+Me 
binary with a period of 12 hours is 2-3R©, which is much 
smaller than the typical size of a red giant (^lOO-lOOOR©). 
This suggests that RX J2130. 6+4710 is a post-common en- 
velope binary (PCEB). In this scenario, the expanding red 
giant star comes into contact with its Roche lobe and be- 
gins to transfer mass to its companion star. This mass trans- 
fer is highly unstable, so a common envelope forms around 
the companion and the core of the red giant. The drag on 
the companion orbiting inside the common envelope leads 
to extensive mass loss and dramatic shrinkage of the orbit 
(Iben & Livio 1993). Common-envelope (CE) evolution af- 
fects the evolution of many binary stars and can produce 
some of the most dramatic objects known, including low- 
mass and high-mass X-ray binaries, novae, supernovae and 
AM CVn binaries. However, it is a poorly understood pro- 
cess, so it is useful to study relatively simple PCEBs like 
RX J2130. 6+4710 to see what they can tell us about this 
phenomenon. 

RX J2130. 6+4710 is one of a very small sample of eclips- 
ing white dwarfs with companions that can be directly de- 
tected in the optical spectrum. These binaries offer a rare op- 
portunity to accurately measure the properties of low mass 
stars and white dwarfs. In this paper we report photomet- 
ric and spectroscopic observations of RX J2130. 6+4710. We 



Table 1. Exposure times, T oxp , for observations of 
RXJ2130. 6+4710 obtained with the JKT. N is the number 
of useful images obtained. 



Filter 


N 


Texp (s) 


U 


6 


120 


B 


1002 


6.5-60 


V 


536 


30 


R 


6 


20 


1 


581 


30 



have used these to measure the masses and radii of the white 
dwarf and the M-dwarf to an accuracy of about 5 per cent. 
We find that the white dwarf in RXJ2130. 6+4710 is a typ- 
ical carbon-oxygen white dwarf and that the M-dwarf has a 
normal radius for its mass. 



2 OBSERVATIONS & REDUCTIONS 
2.1 Photometry 

We first observed RXJ2130. 6+4710 using the lm Jacobus 
Kapteyn telescope (JKT) on the Island of La Palma. We ob- 
tained 478 images with an I band filter and 520 images with 
a V-band filter using a TEK charge coupled device (CCD) 
during the interval 1998 Aug 5-9. The exposure times were 
1.5-3 s for the I band images and 1.5-5 s for the V-band 
images. The image scale on the CCD was 0.33 arcsec per 
pixel. Our intention was to confirm the orbital period of 
the binary from the sinusoidal reflection effect due to the 
irradiation of one side of the M-dwarf by the white dwarf. 
We were delighted to find that the V-band images show an 
eclipse 0.6 magnitudes deep lasting about 30 minutes occur- 
ring at about 0030UT 1998 Aug 10. The I band images show 
an eclipse about 0.1 magnitudes deep at the same time. This 
is exactly what would be expected for the eclipse of a white 
dwarf by an M-dwarf with a similar luminosity in the V- 
band. 

We obtained further photometry of RX J2130. 6+4710 
with the JKT during the interval 2000 August 16-22 us- 
ing a SITe CCD, which gives the same image scale as the 
TEK CCD. Images were obtained using Harris UBVRI fil- 
ters. The number of images secured with each filter and ex- 
posure times are given in Table The shorter exposures 
with the B filter were used for observations during eclipses 
of RXJ2130.6+4710. 

Primary eclipses of RX J2130. 6+4710 were observed 
with ULTRACAM and the WHT on the nights of 17 May 
2002, 24 May 2003 and 25 May 2003. Attempts to observe 
the secondary eclipse on 19 May 2002 and 13 November 2003 
were compromised by poor seeing and flaring of the M-dwarf. 
We have subsequently discovered that the secondary eclipse 
is likely to be extremely shallow (see later) so we do not 
discuss these observations further in this paper. 

ULTRACAM is an ultrafast, triple-beam CCD camera. 
The light is split into three wavelength colours («', g' and r' 
or i') by two dichroic beam-splitters and then passes through 
a filter. The detectors are three back-illuminated, thinned, 
Marconi frame-transfer 1024x1024 active area CCD chips 
with a pixel scale of 0.3 arcsec pixel" 1 . ULTRACAM em- 
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Table 2. Apparent magnitudes of RX J2130.6+4710. 

U B V R I J H K s 



14.8 14.9 14.3 13.70 12.2 11.21 10.58 10.36 



ploys frame transfer CCDs so that the dead time between 
exposures is negligible (24 milliseconds). Signals from Global 
Position System (GPS) satellites were used to ensure that 
the time-stamp for each image is accurate to 1 millisecond, 
(for further details, see Dhillon & Marsh 2001). 

We reduced the data using normal extraction from fixed 
apertures and also Naylor's (1998) optimal extraction from 
apertures varied in step with the seeing. In the latter case 
we used aperture radii 1.5 times the full-width at half max- 
imum of the image of RX J2130. 6+4710. RXJ2130. 6+4710 
is located only 12 arcsec away from a much brighter GO star 
so the photometry can be badly affected in poor seeing. To 
estimate and correct for the extent of this we placed an aper- 
ture on the sky at the same distance from the bright star as 
RX J2130. 6+4710 and symmetrically located with respect 
to the diffraction spikes from this star. 

The bias level in the JKT images was determined from 
the overscan regions and was subtracted from the image 
before further processing. Images of the twilight sky de- 
void of any bright stars were used to determine flat-field 
corrections by forming the median image of 3-5 twilight 
sky images in each filter, one for each night's data. We 
used optimal photometry (Naylor 1998) to determine instru- 
mental magnitudes of the stars in each frame. We checked 
for variability in the stars other than the target star in 
each frame before calculating differential magnitudes be- 
tween RX J2130. 6+4710 and a comparison star located at 
a = 21 h 30 m 19.7 s , S = +47°10'26" (J2000.0). The mean dif- 
ferential magnitudes in the sense (RX J2130. 6+4710 — com- 
parison) in the B, V and I bands when RX J2130. 6+4710 is 
not in eclipse are —0.74, 0.29 and 0.21, respectively. 

We used an aperture with a radius of 5 pixels to measure 
the instrumental magnitude of RX J2130. 6+4710 and the 
5 standard stars in the field of PG 1633+099 with UBVRI 
magnitudes given by Landolt (1992) . We used these to cal- 
culate the apparent magnitudes for RX J2130. 6+4710 given 
m Tabled The calibrat ion of our photometry is approximate 
because the proximity of HD 204906 to RX J2130. 6+4710 
makes it difficult to measure reliable fluxes. For this reason 
we only quote one decimal place for our apparent magni- 
tude values. Also given in Tableware the infrared apparent 
magnitudes for RX J2130. 6+4710 taken from Hoard et al. 
(2002). 

2.2 Spectroscopy 

We observed RXJ2130. 6+4710 with the 2.5m Isaac Newton 
Telescope (INT), also on the island of La Palma, using the 
IDS spectrograph simultaneously with our JKT observations 
in 2000 August. Most of our observations were obtained us- 
ing a narrow slit (0.9-1 arcsec) and a 1200 line mm -1 grat- 
ing on the 235mm camera with an EEV CCD. The spec- 
tra cover the wavelength region 3850-5500A at a resolution 
of about 1.2A and the dispersion is about 0.48A per pixel. 
A typical spectrum is shown in Fig. [5] The exposure times 



were 1200s and we obtained arc lamp spectra approximately 
once per hour to establish the wavelength scale and mon- 
itor the spectrograph drift. We also observed one eclipse 
of RX J2130. 6+4710 with the same instrument but with a 
300 line mm -1 grating and a 4 arcsec wide slit orientated to 
include the same comparison star used for the JKT pho- 
tometry. The exposure times were 10 seconds and there was 
a 4 second delay between exposures to read out the CCD. 
The spectra cover the wavelength range 4215-6795 A at a 
dispersion of 3.7A per pixel. 

Extraction of the narrow-slit spectra from the images 
was performed automatically using optimal extraction to 
maximize the signal-to-noise of the resulting spectra (Marsh 
1989). The arcs associated with each stellar spectrum were 
extracted using the same weighting determined for the stel- 
lar image to avoid possible systematic errors due to the tilt 
of the spectra on the detector. The wavelength scale was 
determined from a fourth-order polynomial fit to measured 
arc line positions. The standard deviation of the fit to the 
18 arc lines was typically 1/20 of a pixel. The wavelength 
scale for an individual spectrum was determined by inter- 
polation to the time of mid-exposure from the fits to arcs 
taken before and after the spectrum to account for the small 
amount of drift in the wavelength scale (< O.lA) due to 
flexure of the instrument between arc spectra. We used a 
spectrum of BD+33 2642 obtained during the same observ- 
ing run and the tabulated fluxes of Oke (1990) to determine 
the calibration of counts-to-flux as a function of wavelength 
in our spectra. No correction for slit losses was attempted 
but we did correct the spectra for atmospheric extinction. 
Statistical errors on every data point calculated from pho- 
ton statistics are rigorously propagated through every stage 
of the data reduction. 

To extract the spectrophotometry taken with a wide 
slit we summed the measured counts over 6 spatial pixels 
centered on the spectrum at every wavelength step for both 
RX J2130. 6+4710 and the comparison star after subtracting 
the sky contribution estimated from a linear fit to the regions 
either side of the spectra. The wavelength calibration was 
determined from a single arc spectrum taken at the start of 
the night and is accurate to about lA. 

2.3 Timing of INT and JKT data 

We warn the reader here that we have not been able to 
establish whether the time-stamps for the data obtained 
with the INT and JKT are reliable. It is probable that the 
time-stamps for these data have an accuracy of better than 
1 second. However, the observatory staff have noted that 
time-stamps on data such as these can be in error by a few 
seconds. 



3 ANALYSIS 

3.1 The ephemeris 

We measured the time of mid-eclipse from our photometry 
using a least-squares fit of a very simple lightcurve model. 
The model is based on the eclipse of one circular disc with 
uniform surface brightness by another similar disc. The pa- 
rameters of the model are the radii of the discs relative to 
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Table 3. Times of mid-eclipse for RX J2130.6+4710. 
BJD (mid-eclipse) Cycle Source 



-1940 -1930 



-1000 
Cycle no. 



Figure 1. Residuals from a linear ephemeris for our measured 
times of mid-eclipse. 



their maximum separation, r± and T2\ the inclination of the 
orbit, i; the orbital period P and time of primary eclipse 
To; the luminosity ratio, l\; the observed intensity outside of 
eclipse, Io and the mass ratio q. We also include a parameter 
to allow for a linear variation of the observed intensity with 
time. As we were only interested in the time of mid-eclipse, 
we fixed the parameters q = 1, i = 90°, P = 0.52104 d 
and found the optimum values of the other parameters us- 
ing the Levenberg-Marquardt method (Press et al. 1992) 
We fitted data covering each eclipse we had observed sep- 
arately. The quality of the fits judged by-eye and from the 
value of x 2 was very good (Fig. [7| . The results are shown 
in Table [3] Corrections from UTC to Barycentric Julian 
date (BJD) were calculated using routines from SLALIB ver- 
sion 2.4-5 (Wallace 2000) The times of mid-eclipse labelled 
"JKT" are from B-band photometry taken with the JKT 
with uncertainties derived from the covariance matrix of the 
least-squares fit. The lightcurve for the time of mid-eclipse 
labelled "INT" was produced by summing the total counts 
over the wavelength 4215-6796A in our low resolution spec- 
tra of RX J2130. 6+4710 and the comparison then taking the 
ratio of these time-series. The uncertainty was also taken 
from the covariance matrix of the least-squares fit. The times 
of minimum labelled "Ultracam" are the weighted averages 
of the values derived by fitting the u , g' and i' lightcurves 
observed with Ultracam. There is an obvious flare during 
the eclipse observed on 2002 May 17 so we have excluded 
data in the range JD= 2452412.632 to 2452412.635 from the 
least-square fit. The uncertainties for these values are taken 
from the standard deviation of the three measurements. The 
measurements from the three lightcurves are consistent to 
the individual uncertainties taken from the covariance ma- 
trix of the least-squares fit. We re-iterate our warning here 
that timing of the data taken with the INT and JKT pre- 
sented here may be in error by a few seconds, so these data 
should not be used to study the long-term period changes. 

A least-squares fit to these data gives the following lin- 
ear ephemeris: 



BJD (Mid - eclipse) 



= (2 452 785.681876 ± 0.000 002) 
+(0.521 035 625 ± 0.000 000 003)£. 



2451775.393803 ± 0.000018 


-1939 


JKT a 


2451776.435891 ± 0.000018 


-1937 


JKT™ 


2451777.477936 ± 0.000005 


-1935 


INT a 


2452412.620375 ± 0.000005 


-716 


Ultracam 


2452784.639806 ± 0.000003 


-2 


Ultracam 


2452785.681873 ± 0.000003 





Ultracam 



All phases quoted in this paper have been calculated 
with this ephemeris. 



"These data are not suitable for long-term period studies. 

3.2 Spectral Type of the M-dwarf 

One of our blue narrow-slit spectra was taken during a total 
eclipse of RX J2130. 6+4710 so it only contains light from 
the M-dwarf. This spectrum is compared to spectra of K- 
and M-dwarfs taken with the same instrument and setup in 
Figure |21 Also shown in Figure |2] is the average of our wide- 
slit spectra taken during eclipse compared to spectra of M- 
dwarfs observed with the same instrument and setup. Spec- 
tral types for the M-dwarfs were taken from Reid, Hawley & 
Gizis (1995) or Hawley, Gizis & Reed (1996). We estimate 
that the spectral type of the M-dwarf in RX J2130. 6+4710 
is M3.5 or M4, based principally on the strength of the TiO 
bands at 5448A 5847A and 6158A and the CaOH band at 
5500-60A(Jaschek & Jaschek 1987). The spectra we have 
used for classification do not cover the TiO bands further 
to the red which are more commonly used for spectral clas- 
sification and we have used only a few stars which are not 
standards for the comparison so the uncertainty in the clas- 
sification is likely to be at least one sub-type. 



3.3 Spectroscopic orbit of the M-dwarf 

We used cross-correlation to measure the radial velocity of 
the M-dwarf from the INT spectra taken in August 2000. 
We first removed the contribution of the broad wings of the 
H/3 and H7 lines from the white dwarf by using two Gaus- 
sian profiles to model each line. The sum of these broad 
Gaussian profiles and a low order polynomial was fit to 
each spectrum by least squares and then subtracted from 
the spectrum prior to performing the cross-correlation. 

We used spectra of four different M-dwarfs taken with 
the same instrumental setup as templates for the cross- 
correlation after normalization using a low-order polyno- 
mial fit by least-squares. The spectra were re-binned onto 
a uniform wavelength grid of 1835 pixels of 32km s^ 1 each 
centered on 4600 A. The cross-correlation function was cal- 
culated excluding 10A around the two Balmer lines, which 
are affected by the emission lines from the cool star and 
the core of the spectral lines from the white dwarf. The ra- 
dial velocity and its uncertainty were measured from the 
cross-correlation function using the method of Schneider & 
Young (1980) which, briefly, determines the centroid of the 
CCF weighted in this case by a Gaussian function with a 
full- width at half-maximum (FWHM) of 200km s" 1 . 

We used a least-squares fit to the measured radial veloc- 
ities, V r , of the function V r = Jm + Km sin(27r[T — To]/P) to 
measure the semi-amplitude of the spectroscopic orbit for 
the M-star, Km- The values of To and the orbital period, 
P, were fixed using the ephemeris described in section |3~T1 
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Figure 2. Left panel From bottom-to-top: spectrum of RX J2130. 6+4710 taken at mid-eclipse, LFT 1580 (M0.5), LFT 1751 (Ml. 5), 
LFT 1466 (M2.5), LFT 93 (M3), mean of 4 spectra of RX J2130. 6+4710 taken near phase 0.75. The spectra have been normalized at 
5100A and offset by 1 for clarity. Right panel From bottom-to-top: mean spectrum of RX J2130. 6+4710 taken at mid-eclipse, LFT 1278 
(Ml), LFT 1436 (Ml. 5), LFT 1580 (M2), LFT 1751 (M2.5), LFT 1466 (M3), LFT 98 (M3.5), LFT 1431 (M4). 



Table 4. Results of least-squares fits of a sine wave plus a distortion to the measured radial velocities of the M-dwarf in RX J2130. 6+4710. 



Template 


Vtcmp 


7M 


K M 


A 


<t>0 


P 


<7"sys 


c 


(Sp. type) 


(kins" 1 ) 


(kins" 1 ) 


(kms" 1 ) 


(kins" 1 ) 






(kins" 1 ) 


(kms- 1 ) 


LFT 1580 


-17.161 


15.6 


135.7 


2.9 


0.18 


0.42 


3.4 


4.5 


(M0.5) 




±0.7 


±0.7 


± 0.8 


±0.04 


±0.07 






LFT 1751 


-27.317 


15.5 


136.4 


2.9 


0.20 


0.40 


3.8 


4.7 


(Ml. 5) 




±0.7 


±0.7 


± 0.9 


±0.03 


±0.06 






LFT 1466 


35.880 


16.1 


136.7 


2.8 


0.23 


0.37 


3.9 


4.8 


(M2.5) 




±0.7 


±0.7 


± 0.9 


±0.04 


±0.07 






LFT 93 


1.500 


13.9 


136.7 


2.8 


0.21 


0.38 


4.0 


4.8 


(M3) 




±0.7 


±0.7 


± 0.9 


±0.04 


±0.06 







From inspection of the residuals from these fits it was clear 
that there was some distortion to the radial velocities with 
an amplitude of a few kms -1 . The residuals tended to be 
positive over the phases [T — To]/P =0.2-0.4 and negative 
over the phases 0.4-0.6. We suspect that this distortion is 
due to a dark spot on the M-dwarf. We do not believe this 
distortion is a direct consequence of the irradiation of the M- 
dwarf by the white dwarf because such a distortion would be 
symmetrical around phase 0.5. To model this distortion we 



add the term Asin([(/> — <j>o]/p) to data in the phase range 
0o < 4> < <t>o + P, where (j> = 2n[T - T ]/P. There may 
be a similar distortion near phase 0.9, but the data near 
this phase range is rather sparse so we have decided not 
to use a more complicated function to fit these data. The 
physical origin of this distortion is a dark spot crossing the 
observer's meridian on the visible disc of the M-dwarf at 
phase m id. The velocity of the spot will vary as sin(^>— </> m id) 
and, since the inclination of the system is approximately 90° 
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Figure 3. Radial velocities for the M-star in RXJ 2130.3+4709 
measured by cross-correlation using the spectrum of LFT 1536. 
The upper panel shows the least-squares fit to the velocities de- 
scribed in the text (sine wave + distortion). The lower panel 
shows the residuals from the sine wave fit and the function fit by 
least-squares which is used to model the distortion. 



its visibility will vary approximately as cos(</> — m id). The 
combination of these effects yields a distortion of the form 
sin(0 — m id) cos(0 — (fimid) = sin[2(0 — m id)]- The effect of 
limb darkening will be to reduce the timescale of the varia- 
tion whereas the finite size of the spot tends to increase the 
timescale of the distortion, so that the velocity distortion is 
approximately Bsin[(<f> — 4>mid)/p] with p ~ 0.5. Note that 
</>mid = 4>a + p/2 so A = —B. 

The value of \ 2 f° r these least-squares fits was still 
rather large. We therefore added an "external error", a syB , in 
quadrature to the standard errors of the radial velocity mea- 
surements to account for additional sources of uncertainty. 
These might be instrumental, e.g., motion of the star in the 
slit, or intrinsic to the M-dwarf, e.g., weak emission lines, 
or spectral features associated with chromospheric activity 
such as flares or dark spots in the photosphere. The value 
of CTsys was chosen to obtain a reduced x 2 value of 1. 

The results of the least-squares fit of the sine wave + 
distortion are given in Table 01 Results are given for each of 
the templates used for cross-correlation as listed in column 
1. The value of 7m has been corrected for the measured 
radial velocity of the template, Kemp, taken from Nidever 
et al. (2002). The spectral type of the template taken from 
Reid, Hawley & Gizis (1995) or Hawley, Gizis & Reed (1996) 
is also given. The standard deviation of the residuals, a is 
given in the final column. An example of the measured radial 
velocities and the function fit by least-squares is shown in 
Fig. El 

We adopt the weighted mean values from Table 0] of 
7m = 15.3 ± 0.7km s" 1 and Km = 136.4 ± 0.8km s" 1 . The 
uncertainties are the standard deviations of the four mea- 
surements with a small additional uncertainty in Km to al- 
low for the partial correlation of Km with A. These values 
are both about 1.5km s" 1 lower than that obtained by fitting 
a simple sine wave to the data. 



3.4 Spectroscopic orbit of the white dwarf 

The measurement of the radial velocities for the white dwarf 
in RX J2130. 6+4710 was difficult because there are only 
three spectral lines available to us which show a sufficiently 
sharp core to give reliable measurements. These spectral re- 
gions all show many sharp absorption lines and a strong, 
variable emission line from the M-dwarf. Fortunately, we 
have a spectrum of the M-dwarf alone obtained during an 
eclipse of the white dwarf. We are able to remove much of 
the contribution of the M-dwarf by subtracting this spec- 
trum from the others. To do this, we first shift the wave- 
length scale of all the spectra to remove the orbital motion 
of the M-dwarf according to the orbit derived in section |3~31 
The spectra were normalized by a constant value determined 
from the mean flux in the line wings. We interpolated the 
spectra onto three wavelength grids centered on each of the 
Balmer lines of 120 pixels of 32km s _1 each. We then formed 
the average of spectra in groups taken at the same orbital 
phase to within 0.02 phase units. 

We then subtracted from these 47 phase-binned spectra 
the spectrum obtained in eclipse scaled by some factor so as 
to optimize the removal of the sharp absorption features. 
This was done independently for each Balmer line as the 
optimum factor varies strongly with wavelength. For the H<5 
line we found that the spectrum obtained in eclipse was 
too noisy, so we used the combined spectrum of the M3 
and M3.5 dwarfs LFT 1431 and LFT 1432 instead, which 
looked very similar but was much less noisy. The resulting 
spectra show the a strong absorption line near the line centre 
and the broad lines of the white dwarf with a weak core 
whose position varies sinusoidally with phase. In addition, 
some spectra show weak emission lines from the M-dwarf 
away from the line centre (Fig. |1J and the stronger central 
emission line from the M-dwarf which is not entirely removed 
by the subtraction process. 

To measure the radial velocity of the white dwarf from 
these spectra, we created a model line profile for the white 
dwarf from the sum of three Gaussian profiles with the same 
mean but independent widths and depths. We then used 
least-squares fitting to the phase-binned spectra to deter- 
mine the widths and depths of these Gaussians. Other pa- 
rameters in the fit were the coefficients of a linear polynomial 
used to model the continuum and the width, height and po- 
sition of a Gaussian used to model the emission line from 
the M-dwarf. Finding the optimum parameters for the Gaus- 
sians and the best scaling factor for the M-dwarf took much 
experimentation, but the quality of the fits achieved judged 
by the x 2 value and by-eye is good. For the measured radial 
velocities reported here, we fixed the widths and depths of 
the Gaussians used to model the white dwarf spectrum to be 
the same for all the spectra. We then repeated least-squares 
fit to find the radial velocity of the white dwarf and the other 
free parameters of the fit. We also corrected for the shift ap- 
plied to the spectra to correct them for the orbital motion 
of the M-dwarf. The phase-binned spectra around the H/3 
line, the model spectra and the residuals from the fit are 
shown in Fig. El Two phase-binned spectra of H/3 observed 
near quadrature and the model spectra are also plotted in 
Fig- El 

We used a least-squares fit to the measured radial veloc- 
ities, V r , of the function V r = 7wd -JfwD sin(27r[T— Tq]/P) 
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Table 5. Results of least-squares fits of a sine wave to the mea- 
sured radial velocities of the white dwarf in RX J2130. 6+4710. 

Spectral Twd %D X 2 CT 

line (kms -1 ) (kms -1 ) (n=47) (kms -1 ) 



H/3 39.2±2.6 139.2 ± 3.4 43.0 17.9 

H 7 58.4±2.7 133.9 ±3.8 43.1 21.4 

H<5 51.7±2.9 135.5 ±4.1 52.7 22.3 



- J 









4830 4840 4850 4860 4870 4880 4890 



Wavelength (A) 

Figure 5. Examples of multiple Gaussian least-squares fits to the 
H/3 line. Two spectra of RX J2130.6±4710 which have had the M- 
star spectrum subtracted offset by 0.4 units are shown (thin line) 
together with the model spectra (thick line). 

to measure the semi-amplitude of the spectroscopic orbit 
for the white dwarf, Kwd- The values of To and the or- 
bital period, P, were fixed using the ephemeris derived 
in section 13.11 Note that 7wd 7m for three reasons: 
i. the gravitational redshift of light from the white dwarf, 
z = O.635(M/M Q )/( J R/R )kms- 1 « 25km s _1 ; ii. the Stark 
effect is asymmetrical for the higher Balmer lines resulting 
in pressure shifts (Grabowski, Halenka & Madej 1987); in. a 
tilt in the continuum will result in a systematic error in the 
value of twd measured from such broad lines. The result- 
ing fit for the HP line is shown in Fig. HJ Also shown are 
the radial velocities measured for the emission line, which 
are a by-product of the analysis. The values of 7wd and 
-Kwd measured from each Balmer line are given in Table |S] 
The agreement between the values of -Kwd measured from 
the three lines is remarkably good, but we suspect this is a 
fluke. These values are certainly affected by systematic er- 
rors due to inaccurate subtraction of the M-dwarf spectrum, 
variable tilts in the spectrum and weak emission lines from 
the M-dwarf. Therefore, we should not take the error in the 
weighted mean as a measure of the uncertainty. Instead we 
use the standard deviation of the values and adopt the value 
-Kwd = 136.5 ±3.8 kms" 1 . 

3.5 The combined spectroscopic orbit 

From the adopted value of P, Km and Kwd we use Ke- 
pler's Law to derive the minimum masses, mm sin 3 i and 
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Figure 6. Radial velocities of the white dwarf (+) and the emis- 
sion line from the M-dwarf (x) measured from the H/3 line. Sinu- 
soidal least-square fits are shown as solid lines. 



Table 6. The combined spectroscopic orbit of RX J2130.6±4710. 



Parameter 


Value 


Notes 


To 


BJD 2 452 785.681876 


See sectionf3.il 


P 


0.521035 625 d 


See section[3.1l 


7WD 


40 kms -1 


Uncertain 


%D 


136.5 ±3. 8km s -1 




7m 


15.3 ± 0.7km s -1 




Km 


136.4 ± 0.8km s -1 




q = mM/mwD 


1.00 ±0.03 




m W D sin 3 i 


0.548 ± 0.016 M Q 




mm sin 3 i 


0.549 ± 0.022 M Q 




asini 


2.81 ± 0.04 R Q 





mwDsin 3 i, the mass ratio q — rriM/wwD, and projected 
semi-major axis asini, where tum is the mass of the M- 
dwarf, mwD is the mass of the white dwarf and i is the 
orbital inclination. The values derived and their uncertain- 
ties are given in Table |S] We have taken the value of 7wd 
derived from the H/3 line as we believe this to be the least 
affected but systematic errors. 

3.6 The lightcurve 

The lightcurve of RX J2130.6±4710 shows that the white 
dwarf is eclipsed for 27 m 30 s and that the ingress and egress 
phases of the eclipse last for lm50s. The depth of the 
eclipse varies from 3.0 magnitudes in the u band to less 
than 0.1 magnitudes in the I band (Fig. 0. There is no sec- 
ondary eclipse apparent in our data, which is unsurprising 
given that it is expected to be less than 1 milli-magnitude 
deep. It is also noticeable that the amplitude of the reflec- 
tion effect is rather small, i.e., the lightcurves are almost flat 
outside of the eclipse. 

For a given inclination, i, a single eclipse will give an ac- 
curate measurement of the radii of the stars relative to their 
separation, i.e., rwD = Rwo/a and tm = Rm/a, where a 
is the separation of the stars and Rwd, Rm are the radii 
of the white dwarf and the M-dwarf, respectively. Addition- 
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Figure 8. The lightcurves of RX J2130. 6+4710 with the following filters (bottom to top); u' , B, g' , V, r' and I. The lightcurves have 
been normalized to the mean out-of-eclipse values and offset by 1.0 for clarity. The data for the u' , g' and r' filters around primary 
eclipse were obtained on the night 2002 May 17 and clearly show two flares. Data at other phases for the u' and g' bands were taken in 
poor seeing on the night 2003 November 13 and also shows two flares. 
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Figure 7. The lightcurve of RX J2130.6+4710 in the u' band (+) 
and the least-squares fit for the eclipse of uniform circular discs 
assuming i = 80.5° (solid line, barely visible). 

ally, the depth of the eclipse gives an accurate measurement 
of the luminosity ratio, l\ at the effective wavelength of the 
lightcurve A. Equally good fits to the eclipse can be found for 
a wide range of i values. The shape of the M-dwarf is accu- 
rately determined by an equipotential surface in the Roche 
potential, which is slightly non-spherical. This gives rise to 
an ellipsoidal effect with an amplitude of about 2 percent in 
the I band lightcurve. The amplitude of the ellipsoidal ef- 
fect varies strongly with vm and the predicted value of vm 
increases as i decreases, so we have enough information in 
the u' band and I band lightcurves to measure rwD, fM and 
i. 

We first measured the values of rwD and ru as a func- 
tion of i from the u band lightcurve using the simple model 
of uniform circular discs described in section 13.11 We in- 
cluded data obtain on the night 2002 May 24 and 2002 May 
25 in the fit after first normalising the data using a straight 
line fit to the data outside of eclipse. The results are plot- 
ted in Fig. [5] We used the slightly more advanced lightcurve 
model ebop (Nelson & Davis 1972, Popper & Etzel 1981) 
to investigate the effect of limb-darkening on the values of 
rwD and tm derived, we find that this radii change by less 
than 0.1 percent if limb-darkening is included in the model. 

We used a simple model to calculate the semi-amplitude 
of the ellipsoidal effect in the I band lightcurve, AJ. The 
model creates a grid of points over an equipotential sur- 
face in the Roche potential. The observed flux as a func- 
tion of phase is then calculated by numerical integration 
of the apparent surface brightness over the points visible at 
each phase. The apparent surface brightness over the star in- 
cludes the effects of gravity darkening and limb-darkening. 
We calculated the value of AJ as a function of i for linear 
limb-darkening coefficients of 0.63 and 0.72, which covers 
the range of values given in the tabulations of Claret (2000) 
for model stellar atmospheres with T e ff =3500 K. The gravity 
darkening exponent was fixed at its standard value of 0.08. 
The effect of gravity darkening is to make the lightcurve 
slightly fainter at phase than phase 0.5. To calculate A J 
we used the semi-amplitude calculated from the flux at phase 
0.5 because the data at phase are affected by the eclipse 
of the white dwarf, although the difference is very small for 



the amplitudes we are dealing with in this case. We also ig- 
nore the effect of the light from the white dwarf when we 
calculate AJ, so a small correction for this contribution to 
the flux has to be made before comparing the predicted and 
observed values of AJ. The predicted values of AJ as a func- 
tion of i are shown in Fig. |U] 

The measurement A J is complicated by the asymmetry 
of the I band lightcurve, e.g., the minimum of the lightcurve 
does not occur at phase 0.5. We assume that this is due to 
the presence of spots on the surface of the M-dwarf. This is 
not an unreasonable assumption given that this is a rapidly 
rotating M-dwarf which is known to show flares (Fig. 8). 
To model this effect of these spots on the lightcurve, we 
assumed that there are two spots which cause dips in the 
lightcurve described by the functions Aj(l + cos[(f> — <pj])/2, 
j = 1,2, which are restricted to orbital phases <f>j ±7r/2. We 
modelled the ellipsoidal effect with the function — AJ cos(2<^)) 
and also included a component for the reflection effect using 
the function — A re f cos($). We first formed the average of 
the observations in phase bins 0.02 wide and calculated the 
standard deviation of the mean of the data in each bin. We 
then subtracted a constant to remove the contribution of 
the white dwarf to the measured flux and renormalised the 
data by its mean. The value of the constant is determined 
from the mean value of the observed flux during eclipse. 

We then performed a weighted least-squares fit of the 
sum of the functions described above plus a constant to the 
phase-binned data excluding the data in eclipse to find opti- 
mum values of the parameters Ai, A2 , 4>i, (j>2, A J, A le { and 
the constant. The resulting fit and the phase-binned data 
are shown in Fig. 1101 With so many parameters, there will 
clearly be many solutions which give a good fit. In this case 
we are interested in finding the likely range of AJ values 
which are reasonable given the observed lightcurve. To find 
the value of AJ we attempted to fit the I band lightcurve 
with the dips due to the spots centered near phases 0.4 and 
0.9 so as to remove as much of the variability as possible 
using the spots alone. The parameters for this fit are shown 
in Table □ 

The I-band data and the radial velocity data for the M- 
dwarf discussed in Section f3 . 31 were obtained simultaneously 
so the parameters for the dark spots we infer by modelling 
distortions to the lightcurve and radial velocity curve should 
be consistent. This is indeed the case for the phase of the 
distortion of the proposed spot with <j>\ « 0.4 which is also 
responsible for the distortion to the radial velocity curve 
with <f>o « 0.2 since <f>i « <f>o + p/2. The radial velocity data 
near phase 0.8 are too sparse to confirm whether the value 
of (f>2 = 0.83 for the second spot is reasonable. We can, in 
principle, check for consistency between the amplitudes of 
the distortions, but this is a less straightforward test (see 
Discussion below). We can at least point out that the size of 
the distortions to the I-band lightcurves are about 2 percent 
and that the distortion to the radial velocity surve is about 
2 percent of the equatorial rotational velocity of the M-dwarf 
(Slkms- 1 ). 

3.7 Masses and radii of the stars 

From the lower panel of Fig. |U] we find that the observed 
value of AJ = 0.0097 ± 0.0024 corresponds to an inclination 
of i = (81.7 ±1.4)°. From the upper two panels we then 
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Figure 10. The phase-binned I band lightcurve of 
RX J2130. 6+4710 with a model fit with two spot-like com- 
ponents, an ellipsoidal component and a reflection component. 



Table 7. Parameters of the least-squares fit to the I band 
lightcurve. 



Parameter 


Value 


Constant 


1.0103 ±0.00053 


AI 


0.0097 ± 0.0024 


A re f 


0.0010 ± 0.0009 


Ai 


0.0175 ±0.0012 


4>i 


0.369 ± 0.023 


A 2 


0.0143 ± 0.0021 


<t>2 


0.833 ± 0.022 


x 2 


40.86 


N 


46 



75 



85 



lnclination( ) 



find rwn = 0.0048 ± 0.0005 and r M = 0.188 ± 0.018. We 
combine these values with the combined spectroscopic orbit 
in Table||j]to derive the masses and radii for the stars given 
in Table |H] The positions of the stars in the mass-radius 
plane is compared to other white dwarfs and M-dwarfs in 

Fig. mi 



3.8 The effective temperature of the white dwarf 

We compared the spectra of RX J2130.6±4710 to a grid of 
synthetic spectra calculated from pure hydrogen model at- 
mospheres to estimate the effective temperature, T e g , of the 
white dwarf. 



Figure 9. Upper panel: The optimum value of rwD as a function 
of inclination from the least-squares fit to the u' band lightcurve; 
middle panel: optimum value of rji; lower panel: the predicted 
semi-amplitude of the I band lightcurve, AI for two values of 
the limb darkening coefficient. The conversion of the observed 
value of AI to values of rwD and tm and their uncertainties are 
illustrated by dashed and dotted lines, respectively. 



Table 8. The masses, radii and effective temperatures of the 
stars in RX J2130. 6+4710. 



Parameter 



White dwarf 



M-dwarf 



Mass (M ) 
Radius (Rq) 
T cB (K) 



0.554 + 0.017 
0.0137 ±0.0014 
18 000 ± 1000 



0.555 ±0.023 
0.534 ±0.053 
3200±100 a 



a Based on our estimate of the spectral type and the calibration 
of Leggett (1992). 
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Figure 11. Upper panel: The white dwarf component of 
RX J2130. 6+4710 in the mass-radius plane compared to other 
white dwarfs as labelled and the models of Benvenuto & Al- 
thaus (1999) for effective temperatures of 7 500K, 17000 K 
and 24 790 K. Lower panel: The M-dwarf component of 
RX J2130. 6+4710 in the mass-radius plane (filled circle) com- 
pared to other M-dwarfs (open circles) taken from Clemens et al. 
(1998). 



These have been generated with the latest versions 
of the plane-parallel, hydrostatic, non-local thermodynamic 
equilibrium (non-LTE) atmosphere and spectral synthesis 
codes TLUSTY (version 200; Hubeny 1988, Hubeny & Lanz 
2003) and SYNSPEC (version 48; Hubeny et al. 1985). All 
calculations include a full treatment of line blanketing and 
use a state-of-the-art model hydrogen atom incorporating 
the 8 lowest energy levels and one super-level extending from 
n=9 to n=80. During the calculation of the model struc- 
ture, lines of the Lyman and Balmer series are treated by 
means of an approximate Stark profile. However, as we use 
the Balmer lines to determine the effective temperature and 
surface gravity of the white dwarf, it is important that we 
employ the best available line broadening data. Therefore, 
during the spectral synthesis step, detailed profiles for the 
HI lines are calculated from the Stark broadening tables of 



Lemke (1997). We calculated synthetic spectra for model at- 
mospheres over a range 15 000 K < T e s < 65 000 K in steps 
of 5000 K and surface gravities in the range 6.5 < log g < 8.5 
in steps of 0.5 (logg in cgs units.). 

In order to compare the spectra of RX J2130. 6+4710 to 
the model spectra, we first had to remove the contribution of 
the M-dwarf from the spectra. The spectra were normalized 
by a constant value determined from the mean flux around 
4600 A. We interpolated the spectra onto a uniform wave- 
length grid of 2760 pixels over the wavelength range 3780A - 
5105A. We formed the average of spectra in groups taken at 
the same orbital phase to within 0.04 phase units. The phase 
ranges were chosen such that spectrum of RX J2130. 6+4710 
taken at mid-eclipse was in a group by itself. We calculated 
a scaling factor to apply to the mid-eclipse spectrum prior 
to subtracting it from each of the 24 phase-binned spectra 
by measuring the ratio of fluxes either side of the band head 
near 4954A and the depth of the absorption feature near 
4226A. The scaling factor calculated was typically 0.25 at 
4226A and 0.36 at 4954A. We multiplied the mid-eclipse 
spectrum by a linear function of wavelength fitted to these 
values prior to shifting it to account for the radial velocity 
then subtracting it from the phase-binned spectra. 

We convolved the synthetic spectra with a Gaussian 
function to account for the resolution of the observed spectra 
and then interpolated the synthetic spectra onto the same 
wavelength grid at the phase-binned spectra. We were then 
able to compare the synthetic and observed spectra by calcu- 
lating the \ 2 statistic for each combination of observed and 
synthetic spectrum. We did this separately for each Balmer 
line after applying a linear normalization calculated from the 
flux in the line wings to all the spectra. We excluded from 
the calculation of x 2 a region ±100kms _1 wide centered on 
the Balmer emission lines from the M-dwarf, the Ca II H and 
K emission lines and an emission line near 3905. 5 A which 
may be due to Si I. 

For 20 of the phase-binned spectra the synthetic spec- 
trum which gave the lowest value of x 2 was calculated 
from model atmosphere with T e g = 20 000 if, logg = 8.0. 
There is good agreement between the value of log g derived 
from fitting these spectra and the value calculated from the 
measured mass and radius, i.e., log<? = 7.93lg Qg. There- 
fore, we created synthetic spectra corresponding to a range 
15 000 K < T efl < 20 000 K in steps of 1000 K and sur- 
face gravities of logg = 7.93, 7.84 and 8.00 by interpola- 
tion over the synthetic spectra calculated from the model 
atmospheres. We calculated the x 2 statistic for these inter- 
polated spectra as before. We found that 17 of the spec- 
tra were best fit by model spectra with T e g = 18 000K, 
log 5 = 7.84 or T cff = 17 000K, log 3 = 7.99. These model 
spectra and a typical phase binned spectrum are shown in 
Fig. 1121 The remaining 7 spectra were best fit by spectra 
within one grid-step of these two models. The quality of fits 
for spectra taken at phases from to 0.5 judged by-eye and 
from the value of \ 2 ls very good. The fits are not so good 
in the other half of the orbit if judged by the value of \ 2 
(reduced x 2 values « 1.5). There is no reason for this which 
is obvious from an inspection of the fits, but it is proba- 
bly due to a combination of inaccurate subtraction of the 
M-dwarf spectrum and inaccuracies in the flux calibration. 
The spectra in the phases range to 0.5 tend to be fit best 
by spectra with T c s = 18 000K. Therefore, we adopt an ef- 
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fective temperature of T e g = 18 000K for the white dwarf in 
RX J2130. 6+4710 and estimate that the uncertainty in this 
estimate is about 1000 K. 



4 DISCUSSION 

We can estimate the distance and reddening to 
RX J2130. 6+4710 from its effective temperature and 
our B and V photometry. From the tabulations of Berg- 
eron, Wesemael & Beauchamp we estimate that the 
white dwarf in RX J2130. 6+4710 has an intrinsic color 
(B— V)o = 0.00 ± 0.02 and an absolute V magnitude, 
Mv = 11.0 + 0.1. From the depths of the eclipses in the 
B and V bands combined with the photometry in Table |2] 
we can estimate apparent magnitudes for the white dwarf 
of V=15.05±0.1, (B-V)=0.2±0.14. From these values 
we estimate E(B-V)=0.2±0.15, (m-M)=4.65±0.5 and a 
distance to RX J2130.6+4710of 85+20 parsecs. 

The distance to HD 204906 is known to be 82+5 parsecs 
from its Hipparcos parallax (Perryman 1987). It would be 
fortunate if RX J2130. 6+4710 could be shown to be physi- 
cally related to HD 204906 because we could then measure 
accurate absolute magnitudes for the stars in this binary. 
Unfortunately, the measured radial velocity for HD204906 
is 34.8+1.4 kms" 1 (Sandage & Fouts 1987), which is very 
different from the value of 7a/ = 15.3 + 0.7 we have mea- 
sured. This suggests that HD 204906 and RX J2130. 6+4710 
are not physically related. One possibility that should be 
investigated is whether HD 204906 is itself a binary star. 

The M-dwarf in RXJ2130. 6+4710 shows all the charac- 
teristics of a magnetically active star, i.e., soft X-ray emis- 
sion, Balmer emission lines, flares and spots. We have at- 
tempted some crude modeling of the distortions to the ra- 
dial velocity and lightcurves to account for these distortions. 
In reality, these spots are likely to be a complex pattern of 
dark regions with variations in the local temperature and 
the emergent spectrum. This pattern of dark regions can 
approximated using models with circular spots, each with 
its own position (latitude and longitude), physical size and 
temperature. Modeling the two spots we have used to model 
the I-band lightcurve would therefore require eight free pa- 
rameters plus some model for how the emergent spectrum 
depends on local temperature and viewing angle. This is 
not warranted given the quality of our data. Further obser- 
vations of RXJ2130. 6+4710 with improved phase coverage 
and resolution could be used to study the pattern of active 
regions on the surface of the M-dwarf and the evolution of 
this pattern using techniques such as Doppler imaging (Rice 
2002) or monitoring of the lightcurve. Such studies would 
benefit from the independent determinations of the rota- 
tional period of the star (assuming it is the same as the or- 
bital period) the inclination of the star and the radius of the 
star presented here. This would, for example, allow a com- 
parison of the rotational period derived from a Fourier anal- 
ysis of the lightcurve excluding eclipses to the orbital period 
in order to look for asynchronous rotation in the M-dwarf. 
Monitoring of the eclipse times will then allow a further test 
of the mechanism proposed by Applegate (1992) to explain 
period changes in binaries similar to RX J2130. 6+4710 such 
as V471 Tau. Indeed, the orbital periods of V471 Tau and 
RX J2130. 6+4710 are very similar, so comparison of the re- 



sults for these two stars will provide observational evidence 
for the way in which Applegate's mechanism depends on the 
mass of the magnetically active star. 

In principle, a measurement of the depth of the sec- 
ondary eclipse in RX J2130. 6+4710 would allow a much 
more accurate determination of the masses and radii of 
the stars. Unfortunately, the predicted secondary eclipse 
depth is 7 milli- magnitudes or less. This will be very dif- 
ficult to measure given the proximity of HD 204906 to 
RX J2130. 6+4710 and the distortions to the lightcurve due 
to spots and flares. A more straightforward way to improve 
the accuracy of the masses and radii will be to obtain im- 
proved spectroscopy, both in terms of resolution and signal- 
to-noise, and to obtain a lightcurve at near infrared wave- 
lengths where the measurment of the ellipsoidal effect will 
be less affacted by star spots. High resolution spectroscopy 
at wavelengths of 3800 - 4000A may reveal sharp metal lines 
due to accretion of metals onto the surface of the white dwarf 
from the stellar wind of the M-dwarf (Zuckermann et al., 
2002). These lines, if present, would allow the measurement 
of accurate radial velocites for the white dwarf and thus an 
accurate measurement of its gravitational redshift. 

The spectral type of the M-dwarf in RX J2130. 6+4710 
is estimated to be M3.5Ve or M4Ve and its mass is 0.555 ± 
0.023Mq. This combination of mass and spectral type is 
not consistent with the calibaration of Leggett (1992) for 
M-stars in the galactic disk. A typical M3.5V in the galactic 
disk has a mass of about 0.2 - 0.3 Mq. An M-type star in 
the galactic disk with a mass of 0.55Mq usually has a spec- 
tral type in the range M0V-M2V. This discrepancy may 
be due to an inaccurate estimate of the spectral type for 
the M-dwarf caused by using a limited spectral range with 
few features sensitive to spectral type and comparing this 
spectrum to stars which are not standard stars for the spec- 
tral class. A much more reliable estimate of the spectral 
type of the M-dwarf in RX J2130. 6+4710 should be made 
using a spectrum taking during the eclipse covering the TiO 
bands in the regions 6000A-9000A. By contrast, the ob- 
served value of (J — H) = 0.63 is negligibly affected by the 
presence of the white dwarf and is typical for M-type stars 
in the galactic disk with a mass of O.55M0. 

We have investigated the past and future evolu- 
tion of RX J2130. 6+4710 using the analysis described 
by Schreiber & Gansicke (2003). The cooling age, t coo i 
of RX J2130. 6+4710 derived by interpolating the cooling 
tracks of Wood (1995) is log(t coo i/y) = 7.85. This is the 
time since RX J2130. 6+4710 emerged from the common- 
envelope phase, at which time the orbital period is estimated 
to have been 0.53- 0.6 d depending on the prescription used 
to model the angular momentum loss due to a magentic 
stellar wind in this interval. The continued loss of angular 
momentum will shrink the Roche lobe of the M-dwarf to the 
point where mass transfer will start from the M-dwarf to the 
white dwarf through the inner Lagrangian point. This will 
occur at an orbital period P s( j = 0.185 d. Our estimate of the 
time before Roche lobe overflow occurs, t B d, also depends on 
the prescription for angular momentum loss used and varies 
from log(t st j/i/) = 8.1 to log(t S d/y) = 9.6. In either case, the 
time taken is less than a Hubble time so RX J2130. 6+4710 
can properly be described as a pre-cataclysmic variable star 
(pre-CV) in the sense defined by Schreiber & Gansicke. The 
value of P sd for RXJ2130. 6+4710 is well above the period 
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gap for CVs and is second only to V471 Tau among the sys- 
tems compiled by Schreiber & Gansicke. This implies that 
RX J2130. 6+4710 is a member of the group of the rare pro- 
genitors of long orbital period CVs. 
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